Electron paramagnetic resonance study of ErSc2N@C8o 
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We present an electron paramagnetic resonance (EPR) study of ErSc2N@C8o fullerene in which 
there are two Er^^ sites corresponding to two different configurations of the ErSc2N cluster inside 
the Cso cage. For each configuration, the EPR spectrum is characterized by a strong anisotropy 
of the ^'-factors {gx,y = 2.9, Qz = 13.0 and gx,y = 5.3, Qz = 10.9). Illumination within the cage 
absorption range (A < 600 nm) induces a rearrangement of the ErSc2N cluster inside the cage. We 
follow the temporal dependence of this rearrangement phenomenologically under various conditions. 



INTRODUCTION 

Rare-earth (RE) doped fullerenes have been widely 
studied because of their interesting structural, magnetic 
and optical properties pji3|. Many studies investigate 
optical properties and are devoted to the determina- 
tion of the electronic structure of the RE ions inside the 
cage [l]|4]-[6]. However, there is a smaller number of re- 
ports about the ground-state spin characteristics of the 
doped ions, which can be examined by electron paramag- 
netic resonance (EPR) |7i (Sj . The combination of EPR 
and optical spectroscopy offers a powerful approach to 
understanding the electronic structure and local symme- 
tries of the ion. 

Photoluminescence (PL) has previously been reported 
from the Er3+ ion in ErScsN^Cgo around 1.5 /im [3l |9l 
[To] , revealing that the ErSc2N cluster inside the Cgo cage 
occupies two dominant spatial configurations, corrobo- 
rating earlier X-ray crystallographic studies [TTj. At low 
temperatures, one of the configurations is favoured over 
the other, and illumination with visible light was found 
to switch the ErSc2N molecular cluster inside the cage 
from one configuration to the other [12j. 

In this article we report a detailed EPR investigation 
of ErSc2N@C8o, including a study on the switching of 
two configurations of the ErSc2N cluster within the cage. 



MATERIALS AND METHODS 

ErSc2N@C8o was purchased from Luna Innovations 
with an approximate purity of 60%. Impurities included 
other erbium doped metallofullerenes. It was further pu- 
rified to >95% by high performance liquid chromatogra- 
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phy (HPLC). Solutions were prepared in toluene (>99%, 
Fisher Scientific) and o-terphenyl (99% Aldrich) with a 
minimum concentration of 3 x 10~^ M (approximately 
10^^ ErSc2N@C8o molecules). The solution was freeze- 
pumped in liquid nitrogen for several cycles to deoxy- 
genate the solvent and was sealed under vacuum in a 
4 mm quartz EPR tube. CW EPR spectra were recorded 
in the temperature range 5 K to 30 K using an X-band 
Bruker EMX Micro EPR spectrometer equipped with a 
rectangular TE102 cavity (quality factor of about 12000) 
with an optical window and an Oxford ESR900 helium 
flow cryostat. The sample was illuminated using a range 
of sources: a Nd:YAG laser (A = 532 nm); a tunable 
infra-red laser (A = 1420 — 1520 nm); and light emitting 
diodes (LEDs) of 10 mW power emitting at 400 nm, 470 
nm and 525 nm. Throughout, the microwave frequency 
was 9.4 GHz. The magnetic field was swept in the range 
to 400 mT, modulated with an amplitude of 0.3 mT at 
a frequency of 100 kHz. 



EPR SPECTRA AND THE HYPERFINE 
STRUCTURE 

Assuming a completely ionic model, the elec- 
tronic configuration of ErSc2N@C8o can be written as 
Er3+(Sc3+)2N3-@C^o [5 • Sc3+, N^- and C^q have no 
unpaired electrons, whereas the Er^+ ion, with 4/^^ elec- 
trons in its valence shell, is presumably the only EPR- 
active species inside ErSc2N@C8o- Erbium has six natu- 
ral isotopes, i^^Er, ^^^^Er, ^^^Er, ^^^Er, ^^^Er and ^^^Er. 
Only ^^^Er, with a natural abundance of 23%, has an 
effective nuclear spin / = 7/2. The other isotopes have 
7 = 0. EPR of rare-earth (RE) ions in a weak crystal- field 
is normally observed only at low temperatures, owing to 
fast spin-lattice relaxation [13]. At 5 K, only the lowest 
doublet is expected to be populated giving rise to the 
resonance spectrum with an effective spin S'eff = 1/2. 
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FIG. 1: (a) X-band CW EPR spectrum of ErScsN^Cgo in a 
frozen toluene solution at 5 K and EasySpin simulated spec- 
trum, (b) Expanded spectra of the experimental and simu- 
lated curves in the low-field region showing the resolved hy- 
perfine splitting. 



TABLE I: Simulated principal values of g and A (in MHz) of 
Er2+ in ErSc2N@C8o for Conf 1 and Conf 2. The and Ay 
values and the components of Q (in MHz) are only reported 
for Er^+ in Conf 1 ; for Er^+ in Conf 2, the effect of these 
values remains unresolved. 





Er2+, Conf 1 




Er3+, Conf 2 


Axis 


g A Q 


9 A 


X 


2.9 (0.3) 380 (50) 


10 (5) 


5.3 (0.3) NA 


y 


2.9 (0.3) 380 (50) 


130 (5) 


5.3 (0.3) NA 


z 


13.0 (0.1) 1380 (10) 


10 (2) 


10.9 (0.1) 1590 (20) 



The EPR spectrum of ErSc2N@C8o in toluene was only 
observed at temperatures below 30 K. The spectrum of 
ErSc2N@C8o in toluene at 5 K is shown in Figure[l] EPR 
peaks, at 51 mT and 220 mT and at 60 mT and 110 mT, 
are attributed to two non-equivalent Er^+ ions inside the 
cage, corresponding to two orientation of the internal 
cluster within the molecule [12 . Additional hyperfine 
structure is also observed around the narrow low-field 
lines. EasySpin software [14 was used to simulate the ex- 
perimental spectrum, with an effective spin-Hamiltonian: 

i^eff = I^BgB.Seii + IJinQuBJ + S^^.A.I + LQ.I (1) 

where the four terms correspond to the electron- and nu- 
clear spin Zeeman interaction, the hyperfine interaction 
and the nuclear-quadrupolar interaction respectively. 

The simulated spectrum, as shown in Figure [l] repro- 
duces the experimental one accurately. The relative in- 
tensities of the EPR peaks are simulated using the nat- 
ural Er isotopic abundance. The simulation parameters 
are given in Table [Tj The A- and Q- tensors are charac- 
terised by one large and two small components, revealing 
the large anisotropy in the system. We extract a set of 



values for Er^+ in each of the two configurations in ErSc2- 
N@C8o- Owing to the narrow linewidth (0.8 mT) of the 
transition at 51 mT for Conf 7, further splitting caused 
by the nuclear-quadrupolar interaction Q was observed. 
However, for the other configurations, Conf 2^ the val- 
ues of Ax and Ay could not be determined accurately 
because of the large linewidths and the low intensities of 
the corresponding peaks of Conf 2. Similarly, the effect 
of the quadrupolar interaction is smeared out over the 
respective hyperfine lines. 

The results of Table |l] seem to suggest an axial sym- 
metry for the tensors. However, from previous X-ray 
diffraction measurements the planar ErSc2N cluster has 
been found to possess a C2v ipi orthorhombic) symme- 
try [11]. In addition, the Er^+ ion sits probably in an 
even lower symmetry site if one considers the neighbour- 
ing carbon atoms of the Cgo cage, which are located 2.5 
to 3 A from the ion [11 . Since the uncertainties obtained 
from our fitting of the A- and Q- tensors are relatively 
large, we then conclude that these tensors have in fact 
no particular symmetry. 

The values of the ratios A^/Qx^ ^y/dy ^z/Qz for 
Conf 1 are not equal. This means that the crystal field 
admixes significant amounts of different {J, mj} states 
to the ground state [13j. The presence of such admix- 
tures has also been shown in the closely related EraN- 
@C8o fullerene from photoluminescence spectroscopy and 
crystal- field modelling [15]. It is interesting that the re- 
sults obtained by Guillot-Noel et al. for Er:YSO crystals 
where the Er sits in Ci symmetry sites, give similar values 
for the hyperfine parameters and thus, significant mixing 
between {J, mj} states [16 . 

As another consequence of the low symmetry of the 
Er^+ in ErSc2N@C8o, the Q tensor axes are not coinci- 
dent with respect to the g and A tensors [17\ The high 
Qy value is about the same order of magnitude as the 
X and y components of the hyperfine tensor, thus mak- 
ing nuclear electric quadrupole transitions observable in 
the EPR spectrum as additional satellite peaks near the 
intense 51 mT peak. 

It is extremely difficult to obtain more structural infor- 
mation from all the tensor values as the EPR spectrum is 
a powder spectrum, further complicated by the presence 
of two configurations whose linewidths are very different. 
Further information could be obtained by aligning these 
molecules, for example with a high magnetic field (some 
magnetic orientation has been shown on Er3N@C8o in 
CS2 with a field of 19 T [15l [18]) and study the EPR 
spectra at various sample rotation angle. 



LINE BROADENING 

At low temperatures, the EPR linewidth 5uj is affected 
by inhomogeneities in the sample such as the g- and 
74-strains, unresolved hyperfine couplings, and also by 
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homogeneous effects which Umit the hnewidths to I/T2, 
where T2 is the spin-spin relaxation time [19 . As the 
temperature rises, the spin-lattice relaxation time (Ti) 
shortens, limiting T2 and leading to broadening of EPR 
peaks. 

Figure |2] shows EPR spectra of ErSc2N@C8o at dif- 
ferent temperatures. The inset shows the evolution of 
linewidths of the peaks at 51 mT and 60 mT as a func- 
tion of 1/T, which fit well to the general expression of 
the form: 



S(jO = 



1] 



Constant 



(2) 



where /c^ is the Boltzmann constant, S is the magnitude 
of the linewidth normalized by the phonon occupation 
number and A^; is the energy of an excited level. This 
dependence reflects that the spin-lattice relaxation is an 
Orbach mechanism, i.e. a two phonon Ti process where 
the relaxation takes place via an excited state [20 . 

The energy A^;, obtained from the fitting parameters 
of Equation ([2| for the two peaks at 51 mT and 60 mT, 
is found to be 36(1) cm~^ and 25(4) cm~^, respectively. 
PL spectroscopy has previously revealed the crystal-field 
energy levels of the Er^+ ion in two different configura- 
tions, corresponding to two different orientations of the 
ErSc2N cluster within the Cgo cage [9 . The energy dif- 
ference between the lowest and second lowest crystal-field 
levels of the ground state ^/i5/2 of Er^+ ion is 37 cm~^ for 
Conf 1 and 28 cm~^ for Conf 2. A^;, deduced from the 
EPR studies, is in good agreement with the energy dif- 
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FIG. 3: Linewidth and g-factor (inset) of Qz feature of 
Conf _/ as a function of microwave frequency. The solid line 
is a fit to a combination of a frequency-independent homo- 
geneous component to the linewidth (12 G), and a g- strain 
component which is linear with frequency (82 mG/GHz). 



ferences obtained in PL. Thus, EPR peaks at 51 mT and 
60 mT are assigned to the optically observed ErSc2N con- 
figurations, Conf 1 and Conf 2^ respectively. 



HIGH-FREQUENCY EPR 

High frequency CW-EPR of ErScsN^Cgo was per- 
formed at the National High Magnetic Field Lab 
(NHMFL), Tallahassee using their home-built 100- 
600 GHz homodyne spectrometer [2T1 . Figure [s] shows 
the linewidth dependence of the narrowest spectral fea- 
ture, the Qz feature of Conf 7, as a function of microwave 
frequency (and thus magnetic field). The data are consis- 
tent with the linewidth at X-band (10 GHz) being homo- 
geneously broadened (on which assumption the analysis 
above relied), while at higher frequency ^-strain across 
the sample causes a linear increase in linewidth. The 
g-strain which can be extracted is small: ~ 0.15% indi- 
cating its value shows a very weak dependence on local 
environment. 



EPR SPECTRA UPON ILLUMINATION 



FIG. 2: EPR spectra of ErSc2N@C8o at different temper- 
atures. The inset shows the temperature dependence of 
linewidths of two most intense peaks at 51 mT and 60 mT ob- 
served in the EPR spectrum at 5 K. The symbols and the solid 
curves indicate the experimental data and the fitted curve, re- 
spectively. The activation energy calculated from fitting the 
evolution of 5uj as a function of temperature support the as- 
signment of these peaks to Conf 1 and Conf 2 respectively. 



Wavelength dependence 

In this section we discuss the influence of optical illu- 
mination on the EPR spectrum of ErSc2N@C8o. For each 
excitation wavelength, the sample was illuminated for 
20 mins at 5 K before the EPR spectrum was recorded. 
Between measurements at different wavelengths, the 
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FIG. 4: (a) Influence of illumination at 400 nm, 470 nm and 
525 nm on the EPR signal of Er^+ ion in ErSc2N@C8o in 
toluene at 5 K. The figures (b) and (c) show changes in 
the peak intensities at 51 mT and 60 mT, and 220 mT and 
110 mT. For a comparison, the EPR spectrum of a non- illu- 
minated sample is also shown. 

sample was heated above 35 K and then cooled down 
to 5 K. This annealing process provided the same initial 
conditions before each new period of illumination. 

Initially, we illuminated the Er^+ ions to excite tran- 
sitions within the crystal-field levels directly, using ex- 
citation wavelengths of 1496 nm and 1499 nm [9] for 
Conf 1 and Conf 2 respectively. These wavelengths are 
too long to excite the cage and had no effect on the EPR 
spectrum. 

We then excited ErSc2N@C8o at different wavelengths 
well within the cage absorption range using 400 nm, 
470 nm and 525 nm LEDs. The EPR spectra of ErScs- 
N@C8o under LED illumination are shown in Figure |4] 
For comparison, the EPR spectrum without illumination 
is also shown on the same graph. The spectra reveal that 
the intensities of the EPR peaks corresponding to the two 
ErSc2N configurations ions vary in opposite directions 
(indicated bt the arrows), allowing the assignment of the 
set of peaks at 51 mT and 220 mT to Conf 1 and the set 
of 60 mT and 110 mT peaks to the other Conf 2 [12]. 
The line shape and spin-Hamiltonian parameters remain 
unchanged on illumination, and no evolution of any new 
spectral features is observed. 

The change in the EPR intensity for each configuration 
indicates that the ErSc2N cluster rearranges within the 
Cgo cage upon illumination. Exciting ErSc2N@C8o with 
visible light appears to preferentially drive the switching 
one configuration to the other. As the direct excitation 
of Er^+ ion has no influence on the relative proportion 
of the two configurations, we can conclude that the con- 
figurational changes are mediated via the cage. Subse- 
quently, we used a green laser operating at 532 nm as an 
excitation source. The resulting change in the EPR spec- 
trum was consistent with the effect of illumination using 
LEDs. Probably owing to the higher excitation power 
(26 mW, compared to 10 mW for the LEDs), illumina- 




FIG. 5: Populations of Conf 1 and Conf 2 from the ESR 
peaks of ErSc2N@C8o in o-terphenyl for various experimental 
conditions. The exictation laser was set to 532 nm and 20 mW 
power, except for one measurement as indicated in the figure. 
The sample was annealed at room temperature between each 
measurement. 

tion with the 532 nm laser induced a 10% larger change 
in the EPR intensities as compared to the illumination 
with the 525 nm LED in the given time span of 20 mins. 

Temperature dependence 

We studied the effect on the EPR spectrum of laser 
illumination at 532 nm at various temperatures, T, be- 
tween 5 K and 15 K. The integrated area of the spectrum 
depends on many factors, including transition moment, 
temperature, microwave power and relaxation times. We 
combine all such terms into a constant of proportionality 
K such that the intensity of a feature attributed to Conf z, 
li = hZiUi^ where is the number of spins in Conf i. 

Assuming that HCi for a particular configuration is inde- 
pendent of the illumination, the change in the integrated 
area under the resonance peak provides a measure of the 
change in the configuration populations. We examined 
the resonance peaks at 51 mT and 60 mT corresponding 
to Conf 1 and Conf 2 respectively. Figure [5] shows the 
change in the populations for Conf 1 and Conf 2 at var- 
ious temperatures with and without illumination. The 
populations were obtained by a fitting procedure assum- 
ing that the total number of spin is constant and that 
only two configurations are present. The latter assump- 
tion has been verified experimentally (see the population 
switching dynamics section). 

Previously, it has been shown that at 5 K a swing in 
configuration population by as much as « 27% Conf ^is 
possible [12 upon optical illumination. In order to study 
the effect on switching of various parameters, such as 
temperature, it is important to control the molecular 
environment carefully. The use of toluene as a solvent 
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produces a metastable glass upon freezing which eas- 
ily cracks and form a polycrystal. In such a polycrys- 
tal, the various orientations of the grains induce strains 
that probably impact on the local environment of the 
fullerene, and hence the Er^+ ions. We therefore used 
o-terphenyl as a solvent as it is known to produce stable 
glasses, and therefore a more reproducible environment 
for the ErSc2N. In this case, the sample needs to be first 
heated to about 330 K and then quickly frozen in liquid 
nitrogen. We observe in figure [5] that the Conf 2 is pre- 
dominant at higher temperatures. Moreover, a weaker 
laser illumination results in less switching of Conf 1 to 
Conf 2. To summarise, the analysis of the populations 
before and after illumination shows that higher temper- 
atures and illumination both favour Conf 2. 



ated. As we know that: 



/i(oo) ni{oo)tzi 
12(00) 712(00)^1:2 



(3) 



and the ratio ^^2/^1 = 3.2, we can deduce the ratio 
ni(oo)/n2(oo). At 5 K and 20 mW, the ratio is: 



ni(oo) 
^2(00) 



: 0.45(2) 



(4) 



This ratio is in agreement with that reported previ- 
ously [12]. We note that this shows the light-induced 
switching observed is not due to sample heating, which 
would reduce the intensity of both peaks, as shown above. 



POPULATION SWITCHING DYNAMICS UPON 
ILLUMINATION 

Evidence for the switching between only two 
configurations ErSc2N cluster 

Here we report the experimental dynamics of the pop- 
ulation switching between two configurations of ErSc2N 
cluster within the Cgo cage upon illumination. 

The magnetic field range was selected to include the 
two sharp peaks at 51 mT and 60 mT. The sample was 
diluted in o-teprhenyl as it formed a better glass than 
toluene. After being kept in dark for about 5 min, it was 
illuminated with a laser at 532 nm and 20 mW power for 
48 min. The changes in EPR peaks at 51 mT and 60 mT 
were recorded simultaneously. 

As the laser is turned on, the intensity of the peak 
at 51 mT, corresponding to Conf 7, decreases and the 
one for the peak at 60 mT, corresponding to Conf 2^ 
increases. 

The integrated area (7^) under the two sharp peaks 
at 51 mT and 60 mT was calculated for each spectrum. 
Figure [6]^a) shows the temporal evolution of the areas 
{li) before and during illumination. 

We assume that the total number of spins ni + n2 = 
— + — is constant throughout illumination. We test this 
assumption and derive the ratio of the intensity factor z^^, 
by plotting {^Ii -\- 12) as function of time for different 
ratios of 1^2/ 1^1 • Through a fitting procedure we find that 
for /^2//^i = 3.2 ± 0.1, the function {^h + ^2), which is 
proportional to the total number of spins in Conf 1 and 
Conf 2^ is constant. This is illustrated in Figure [6]^b) 
for a range of 1^2/1^1 = 2.2,3.2,4.2. This value can be 
approximately understood by considering the different 
total spectral widths {g± to g\\) of the two configurations. 

As t ^ 00, the area under the two peaks corresponding 
to each configuration, /i(oo) and 12(00), reaches a steady 
regime under a constant illumination (see Figure [6]^ a)). 
From Figure [6ja), the ratio /i(oo)//2(oo) can be evalu- 



Eff*ects of temperature and power on the temporal 
evolution 

In order to gain further insight on the switching dy- 
namics, we performed similar time scan measurements 
at a lower power and higher temperatures. The tempo- 
ral evolution of the integrated area under various con- 
ditions for each configuration is shown in Figure [7| We 
confirmed that for most conditons, the ratio K.2/1^1 for 
3 mW is unchanged at 3.2 ±0.1. Only at 10 K this ratio 
is slighlty higher, with K2I — 3.8±0.1. Thus, we can 
directly translate the changes in EPR peak areas to the 
populations dynamics of the configurations. 

The data are well fit to biexponential decay functions 
of time, with a slow decay time (T5) and a fast decay 
time (Tj). We consider the average decay time Tavg = 
{AsTs-\-AfTf), where Ag and Af are the normalised am- 
plitudes of the two exponential decay components. The 
results are presented in Table |n[ Reducing laser power 
from 20 to 3 mW has no significant effect on the ex- 
tracted ratios of ni(0)/n2(0) and ni(oo)/n2(oo), as ex- 
pected from Figure 5, however the average decay time 
Tavg is much slower. 

In contrast, changing the temperature from 5 to 10 K 
involves a change in both the initial and final population 
ratios (as observed above). However it does not affect 
Tavg- Thus, the temporal evolution of the switching pro- 
cess is dominated by the kinetic effect of laser power, 
rather than a thermodynamic effect. 

Simple kinetic models involving different excitation 
and relaxation rates can reproduce the basic form of this 
biexponential behaviour, however a quantitative compar- 
ison would require a detailed understanding of the en- 
ergy transfer process between the carbon cage and the 
ion . Further optical studies on ErSc2N and related 
TNT fullerenes will help to elucidate this transfer process 
and the energy levels involved. 
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FIG. 6: (a) Temporal evolution of the integrated area under the EPR peaks at 51 mT and 60 mT before and during illumination. 
Temperature = 5 K; excitation wavelength A = 532 nm; power = 20 mW. (b) Plot of 1^2/1^1 * A +^2) versus time for 1^2/1^1 =2.2, 
3.2 and 4.2. 
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FIG. 7: (a) Time response of the integrated area under the peak at 51 mT before and during illumination at T = 5 K and 10 K 
with an excitation wavelength A = 532 nm and power of 20 mW and 3 mW. (b) Time response of the integrated area under 
the peak at 60 mT before and under illumination at T = 5 K and 10 K with an excitation wavelength A = 532 nm and a power 
of 20 mW. The laser was turned on after 400s. Smooth curves show biexponential fits with the parameters from the model. 



T(K) 


Power (mW) 


ni(0) 
n2(0) 


ni (00) 
722(00) 


Tavg (s) 


5 


20 


1.33(0.02) 


0.45(0.02) 


302(60) 


5 


3 


1.28(0.02) 


0.50(0.02) 


613(100) 


10 


20 


1.12(0.02) 


0.37(0.02) 


309(60) 



TABLE II: Average decay times deduced from the fitting 
parameters of the bi-exponential temporal evolution of the 
average integrated area as a function of temperature and laser 
power. 



Relaxation after illumination 

When the laser is turned off, the peak intensities evolve 
in the opposite direction in a monoexponential way. This 
is consistent with a direct relaxation between Conf 2 and 
Conf 1. This relaxation is very slow compared to the 
excitation timescales: after 12 hours at 15 K, the con- 
figurations relax only half way to equilibrium. Thus at 
low temperatures, the configuration populations are very 
stable, while they recover to equilibrium very rapidly if 



the sample is heated above 35 K. 



CONCLUSIONS 

We have investigated the behaviour of two dominant 
configurations, Conf 1 and Conf 2^ of the Er^+ ion in 
ErSc2N@C8o using continuous wave EPR. The experi- 
mental spectra can be simulated to extract a spin Hamil- 
tonian with anisotropic g- and A- and Q- tensors for each 
cluster configuration, to the extent allowed by the EPR 
linewidths. The results suggest that the erbium ion sits in 
a low symmetry site, at least for Conf 1. The two config- 
urations are switchable through optical excitation of the 
Cgo cage. The switching dynamics can be described phe- 
nomenologically by biexponential functions of time. The 
average decay time becomes slower at lower laser power 
but does not change with temperature. The recovery to 
equilibrium populations is very slow at low temperatures, 
but the equilibrium populations can be quickly recovered 
by annealing the sample above 35 K. This thermo-optical 
hysteresis makes ErSc2N@C8o a potential candidate for 
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reversible switches and memory elements [22H28] . 
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